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Abstract 

Many networked systems involve multiple modes of transport. Such systems 
are called multimodal, and examples include logistic networks, biomedical 
phenomena, manufacturing process and telecommunication networks. Exist- 
ing techniques for determining optimal paths in multimodal networks have 
either required heuristics or else application-specific constraints to obtain 
tractable problems, removing the multimodal traits of the network during 
analysis. In this paper weighted coloured-edge graphs are introduced to 
model multimodal networks, where colours represent the modes of trans- 
portation. Optimal paths are selected using a partial order that compares 
the weights in each colour, resulting in a Pareto optimal set of shortest paths. 
This approach is shown to be tractable through experimental analyses for 
random and real multimodal networks without the need to apply heuristics 
or constraints. 

Keywords: graph theory, multimodal network, Pareto optimal set, shortest 
path, weighted coloured-edge graph 



1. Introduction 

Extensive scientific literature has been devoted in the last three decades 
to the study of multimodal networks (MMN). During this time, research 
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has main ly focused on practic al ap plications for freight or urban tran sporta 



tion (see IJarzemskiend ( 120071 ) and iMacharis and Bontekoningj ( 12004| ) for ex- 
tensive reviews). As a system in which several means of transport are 
available, a multimodal system is able to emulate a wide spectrum of real 
life phenomena beyond the field of logistics. Areas such as computer net- 
works, biology and manufacturing have begun to utilize multimodal net- 
works f or studying ari d mod el ling situations. Ex a mples can be found in pa - 
pers bv lAbrach et all f[2003h. IChen et all fl2005h. iHeath and Siosonl fl2007h 



Medeiros et all (12000f l 



Kiesmiiller et al.l (120051 ). iNigay and CoutazI (119931 ) 
and I Siosonl (120051 ). 

From a modelling perspective, a range of techniques have been used to 
model MMN. They can be loosely classified into three predominant domains: 
mathematical programming, weighted graphs and multi-weighted graphs. 

1.1. The mathematical programming approach 

These techniques are characterized by making use of linear or non-linear 
formulations for representing a MMN by a set of equations. 

Linear programming techniques are suitable when each decision variable 
i s a li near combination of the problem parameters, iHillier and Lieberman 
( 120091 ). Integer programming and mixed integer programming stand out 
as the most common linear programming techniques used for multimodal 
modellin g. Sample p aper s using linear pro gramming as a modelling tool are 
given bv lMinI (Il99l[ ) and lKim etaP ( 1l999h . 

Non-linear programming is another renowned modelling technique for 
MMN. It is mainly used to build intricate cost functions, and principally 
deals with second order eq uations satisfying con v ex or c oncav e pro perties 



Exarn ples are provided by iKim and KimI (j2006[ ). iHornI (j2003[ ) and IChang 
( 120071 ) which have opted to use non-linear programming as their main mod- 
elling approach. 

In the mathematical programming approach, mode options are visualized 
as decision variable indices, which considerably increases the complexity of 
the problem. Relaxation or cutting plane techniques are commonly used to 
make the problem tractable. Interesting papers tackling general views of 
mathematical programming for intermodal transp ortation (the t ransp orta- 
tion of goods) a nd urban transportation were by IJarzemskiend (120071 ) and 
Nagurneyl ( 119841 ) respectively. 
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1.2. The weighted graph approach 

In this approach, a node typically represents a location, such as a ware- 
house, transportation hub or network router, and an edge represents a trans- 
portation link, such as rail line, a bus or a wireless connection. A variety of 
graphs have been used to study these tra nsport systems, su ch as digraphs, 
multigraphs, hypergraphs and grid graphs. lAyed et al.l ( 120081 ) provides a gen- 
eral classification for MMN models based on weighted graph approaches. In 
particular, the article emphasizes the use of multigraphs, in which there might 
be multiple edges between two nodes, and the use of grids in which a grid 
is overlayed on a planimetric map. Both can result in dense graphs, which 
require edge reduction techniques to make their analysis tractable. In prac- 
tice such reductions rely on enforcing constraints on feasible edges in order 
to build a specifi c path. St u dies making use of such grap hs are provided by 
Foo et all (Il999h . lOiang Lil (l2000l ). iKitamural (Il999h and iFragouli and Dehs 
f l2002h . Hypergraphs are another type of graph used in some articles. In 
graph theory, a hypergraph is a general i zation of a graph, where edges can 
connect any number of vertices, iLawlerl (120011 ). In the multimodal context, 
such graphs have found interesting applications in biology and urban trans- 
po rtation. Sample papers u s ing hypergrap hs to represent MMN are yielde d 
by iHeath and SiosonI (120071 ) , ISiosonI (120051 ) and iLozano and Storchil (|2002[ ) . 

In effect, the weighted graph approach only utilizes mode information dur- 
ing the application of constraints, removing the multimodal traits from the 
network during modelling. The analysis in this approach is very application- 
dependent as it relies on applying application-specific constraints. 



1.3. The multi-weighted graph approach 

This approach has been extensively utilized for the Multicriteria Shortest 
Path Problem (MSPP) which h as be come an f r uitful branch of research since 
the 1980s, see iTarapatal (120071 ) and ISoroushI (120081 ) for a complete review. 



Basically, the approach assigns multiple weights to each edge. In particular, 
the bicriteria shortest path problem assigns two weights to each edge, such 
as cost and time. 

Optimality in the multi-weighted graph approach is commonly estab- 
lished by the use of a partial order relation which results in a Pareto op- 
timal set of paths that are candidates for the sought shortest path. There 
is little literature that directly applies multi-weighted graphs for modelling 
MMN, but the goal of MSPP is essentially the same as for the shortest 
path problem in MMN. Although articles developing MSPP formulations 
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for MMN can be identified, they preferentially use partial orders to op- 
timize route choice decisions mainly based on cost and time, leaving the 
mode options as an outcome of the optimal route. The tractability of the 
MSPP is i nextricably connected with the cardin a lity of the Par eto set ac- 
cording to iMiiller-Hannemann and Weihd (|2006[ ) . iHansenl (jl980( ) showed in 
his paper t hat this cardinality is exponential in the worst case, although 



Louil (Il983[ ) points out that Pareto sets for some graphs with multidimen- 
sional weights have polynomial average case cardinalities. Constraints are 
applied during analysis to make the problem tractable, resulting in a Pareto 
set with manageable cardinality. MMN mod els whose mainstay is a multi^ 
weigh te d graph can be found in pape rs by lAndroutsopoulos and Zografos 
fl2009h . lAifadopoulou etaP fcoOTh and lModesti and Sciomachenl dlQQSh . 

The paper is organized as follows. In Section 2 a fourth approach for mod- 
eling MMN is introduced and compared with the three previous approaches. 
An algorithm is described in Section 3 that is used in the paper to compute 
shortest paths in weighted coloured-edge graphs. In Section 4 the model 
is experimentally studied to demonstrate the tractability of the approach. 
Moreover, the algorithm is applied to a real multimodal network in order 
to assess its practical applicability. Finally, Section 5 provides conclusions 
about this research. 



2. Weighted coloured— edge graph approach 

All the approaches described in Section [T] are heavily application specific 
and do not actually utilize the multimodal nature of a network. In this paper 
an approach to model and analyze multimodal networks is introduced. In 
essence, such approach uses a weighted graph in which edges are endowed 
with two attributes: a positive weight and a discrete variable called colour. 

A weighted coloured-edge graph G = (y,E,co,X) consists of a directed 
graph (V, E) with vertex set V and edge set E, a weight function u: E R"*" 
on edges, and a colour function X: E M on edges, where M is a set of 
colours. Typically M is taken as a finite set with k = \M\. Associated 
to each edge e„„ G E from vertex u to vertex v there is a positive weight 
uj{euv) and a colour X{euv)- For any colour i E M and for any path = 
{exoxi, ex'ixa, 62:2x3, • • ■ , (^xi_-,xi} between two vertices u = Xq and v = Xi, where 
each Xi G V, the path weight C0i{puv) in colour i is defined as Ui{puv) = 
Se,,,^+,ep,A(e,,,^+i)=i^(ex,x,+J- The total path weight is represented as a k- 
tuple {uji{puv), 0Ji{puv), uJkiPuv)), giving the total weight of the path 
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in each colour. 

Note that there is no restriction placed on the number of edges Cuv from 
a vertex to a vertex v. However, in practice for the shortest path problem 
attention can be restricted to weighted coloured-edge graphs for which there 
is at most one edge e„„ in each colour from u to v. 

Let u and v be two given vertices of G and let be the set of all paths 
from u (source) to v (destination) in G. A binary relation between two paths 
Puv and is defined by puv < Puv if and only if Ui{puv) < oJiipL) for all 
i. The relation < is clearly reflexive and transitive and gives a partial order 
on the A;-tuple path weights, but only a preorder on the paths themselves as 
multiple paths might have the same total path weight. 

Let Muv = {Puv e ^uv I ^p'uv ^ "^uv with u{p'^J 7^ u{puv),^i < k such 
that Ui{puv) < ^i{Puv)} be the set of Pareto optimal paths joining vertices 
u and V. This set has an important characteristic: for any p^v G Aiuv, it 
is impossible to determine a path p'^^ from u to f which has smaller weight 
than Puv in some of its k c olours wi t hout at least one of the other weights 
being larger, analogously to iMartinsI (Il984[ ). 

From the above definitions it is apparent that the concept of a weighted 
coloured-edge graph with k colours can equivalently be formulated as a mul- 
tiweighted multigraph where each edge is assigned a fc-tuple of non-negative 
weights {wci, ■ ■ ■ ,Wc^, . . . ,WcJ and exactly one > 0. However, multi- 
weighted graphs are mostly used in multicriteria optimization applications 
where the weight components correspond to quantities to be optimized, such 
as cost and time, so edges typically contribute toward more than just one 
quantity. For this reason multiweighted graphs whose edge weights are zero 
in all but one component have not received attention in the literature. 

Shortest path analysis in the weighted coloured-edge graph approach is 
seen in this paper to typically be tractable without the need to apply any 
application-specific heuristics or constraints, so can be considered a general 
tool for the study of multimodal networks. Application-specific considera- 
tions can still be applied to the resulting set J^uv, or a post-optimal analysis 
undertaken on it. One facet of this model is that it can be directly applied 
to multigraph applications, such as transportation networks where there are 
multiple transportation means between two locations, communication net- 
works where there are multiple links or choice of communication protocols 
between nodes, or epidemic models which have multiple paths of infection. 

However, focusing attention on only the Pareto optimal paths limits the 
approach to shortest path problems where just the summed contribution of 
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each colour is important, and where any measure of optimahty is presumed 
to be an increasing (hnear or non-hnear) function of the summed contribu- 
tion in each colour. For instance, the approach presumes in a transportation 
network that the optimal path (such as least cost, time, or distance) is some 
application-specific increasing function of the total weight in each transporta- 
tion means, or that the user can apply some application-specific criteria to 
select a preferred path from the Pareto set once it has been determined. 

The approach can be adapted for path constraints such as restricting 
the number of hops or the number of mode changes by slightly enhancing 
the algorithm used to determine the Pareto set. For instance, besides using 
colours to represent the different transportation means, an additional colour 
can be used to count the number of edges in a path as the path is being built 
during the analysis and/or to count the number of transfers from one means 
of transportation to another. 

Optimization problems that utilize models similar to weig hted coloured 



edge g raphs have received little attention in the literature. iClimaco et al 
(I2OIOI ) experimentally studied the number of spanning trees in a weighted 
graph whose edges are labelled with a colour. In that work, weight and colour 
are two criteria to be both minimised and the proposed algorithm generates 
a set of non-dominated spanning trees. The comput ation of co l oured paths 
in a weighted coloured-edge graph is investigated in IXu et al.l (120091 ) . The 
main feature of their approach is a graph reduction technique based on a 
priority rule. This rule basically transforms a weighted coloured-edge mul- 
tidigraph into a coloured-vertex digraph by applying algebraic operations 
to the adjacency matrix. Additionally, the authors provide an algorithm 
to identify coloured source-destination paths. Nevertheless, the algorithm is 
not intended for general instances because its input is a unit weighted colored 
multidigraph and only paths not having consecutive edges equally coloured 
are considered. 

This paper investigates the feasibility of the approach as a general tool for 
multimodal networks by determining the cardinality of the Pareto optimal 
set Aiuv for many diverse networks. It shows that the cardinality is typically 
a very low order polynomial function of the size of the network, and demon- 
strates that even dense multimodal graphs with hundreds of thousands of 
edges can be feasibly analyzed using this approach, without the need for any 
reduction techniques. In fact, it is seen that the number of modes k is more 
of a limiting factor of the approach than is the number of vertices or edges 
in the graph. 



6 



3. Algorithm for determining Pareto Optimal Sets 



To experimentally study the feasibility of using weighted coloured-edge 
graphs for multimodal networks an algorithm that determines Muv is re- 
quired. A simple generalization of Dijkstra's algorithm from unimodal net- 
works has been developed for the purposes of this paper, although more ef- 
ficient algorithms might be investigated in the future. The classic Dijkstra's 
algorithm for solving the single-source shortest path problem in unimodal 
networks uses a priority queue Q to store shortest path estimates from a 
fixed source vertex s to each vertex v in the network until the shortest path 
to V is determined. Since the weights of any paths Ps^j from s to v are linearly 
ordered there is only at most one shortest path estimate in the queue at a 
time for each vertex v. At the start of each iteration of the algorithm the 
shortest path estimate at the front of the queue is the actual shortest path 
to one of the vertices in the network. 

In a weighted coloured-edge graph Dijkstra's algorithm must be slightly 
generalised to handle weights of paths being partially ordered rather than 
linearly ordered. A priority queue Q can again be used to store shortest path 
estimates with the requirement that if a path psv from s to v has smaller 
weight than another path then it must appear earlier in the queue. Al- 
though the results presented in this paper use such a simple queue instead of 
a more sophisticated non-linear data structure the performance of the algo- 
rithm is seen to be surprisingly good. As in the classical Dijkstra's algorithm 
the weighted coloured-edge version of the algorithm takes as input a network 
G and a source vertex s. It commences at s with the empty path pss and 
relaxes each edge that is incident from the source vertex s, adding the single 
edge paths to the queue. At the front of the queue will be a shortest path 
estimate psv to some vertex v adjacent to s. Since all weights are positive 
in the network psv must have minimal weight amongst paths from s to v 
(although it might not be the only minimal path from s to f in the queue), 
so Psv is added to the set A4sv and removed from the queue. The algorithm 
then relaxes all the edges incident to v, extending the path pgy by each edge 
to a path p'g^ = Psv U {cvu}, adding those extended paths p's^ to the queue 
that have minimal weight amongst paths from s to u, and removing from the 
queue any path p"^ from s to m that has greater weight than p'^^. The algo- 
rithm repeats itself until the queue is empty, producing as output the Pareto 
optimal set Aisv for each vertex v in the network. Figure H] describes the 
pseudocode of the algorithm using the notation developed by 



Gormen et al. 
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Multimodal-Dijkstra(G, s) 

1 [> Initially no Pareto optimal paths known 

2 for each vertex v 

3 do Msv ^ 

4 > Create a queue Q to hold shortest path estimates during processing 

5 g^0 

6 add the empty path from s to s into Q 

7 while g 



8 do remove the path at front of Q that has some end vertex v 

9 add the path psv to A4sv 

10 > Relax the edges incident from v 

11 for each edge e^u incident from f to a vertex u not in 

12 do > Extend the path psv by the edge e^u 

13 ^ Psv U {e„„} 

14 if p'g^ has minimal weight in Q from s to 

15 then add p'^^ to Q 

16 > Remove any paths no longer minimal in Q 

17 for each e Q with u {p'^J > w {p'^J 

18 do remove the path p"^ from Q 



19 return M 



Figure 1: Pseudocode of the algorithm 



(120011) 



The number of relaxation steps is an important indicator of the algo- 
rithm's order, so besides finding Mlsv the experiment discussed in Section H] 
also track the number of paths p'^^ processed by the algorithm. 

As an example of an application of the weighted coloure d-edge graph ap- 



proach , the algorithm is run with a multimodal network from lLozano and Storchi 



(120011 ) starting at source vertex 0. Figure [2] shows the network which has 
21 vertices, 51 edges and 4 different transport choices (bus, metro, private 
and transfer). The algorithm commences with just the empty path poo on the 
queue and relaxes two edges: eoi with weight (bus, metro, private, transfer) = 
(15,0,0,0), and Cqs with weight (0,0,5,0), which are both added to the 
queue. Since the two weights are incomparable, either could be at the front 
of the queue, so the next iteration of the algorithm either adds the path 
Poi = {coi} to Aioi and relaxes the four edges incident to vertex 1 by ex- 
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tending the path poi by each, or else adds the path po3 = {^03} to A^oa 
and relaxes the three edges incident to vertex 3 by extending the path pos by 
each. Continuing in this way the Pareto optimal set Aiov is obtained for each 
vertex v in the network, resulting in 52 Pareto optimal paths from vertex 
to vertex 20 whose weights are listed in Table [1] Depending on the applica- 
tion, constraints or heuristics can then be applied to the 52 paths to select a 
path preferred by the user. Using just a simple priority queue data structure 
the generalised Dijkstra's algorithm can det ermine Aiov for all 21 vert ices 
V within approximately 10ms. The article of iLozano and Storchil (120011 ) in- 
stead uses a weighted graph approach with application-specific constraints 
and a simple cost function which adds the weights in each mode together to 
get a single valued total weight, resulting in the paths numbered 2, 25, 33, 
47 in the table. 

Note that a Pareto set permits a post-optimal analysis to be carried 
out provided that the total cost is presumed to be an increasing function 
of the weight in each mode. For instance, suppose in the example net- 
work that the edge weights represent distance and for simplicity presume 
the cost is one dollar per unit distance for each means of transport. A 
natural optimization question could be how much the unit cost associated 
to a particular mode could be increased or decreased with the current op- 
timal path remaining optimal. As an illustration, path 25 which has the 
edges {eo3, 631, Cig, eg 10, Ciom, Cmis, 61517, e^ie, Cieis, eisig, eig2o} and shown 
in Figure E] has least total cost $47, but from the Pareto set it is easily seen 
that an increase of over 20% in the relative metro costing would make path 
41 with edges {cqs, 631, eig, 6913, 61315, 61517, 617 le, eie is, eis 19, eig 20} a better 
choice, or a 25% increase in bus prices would make path 16 with edges 
{^03, 632, 6210, ^1014, ^1415, 61517, 61716, 616I8, Cisig, 61920} better. 

This example demonstrates that the Multimodal Dijkstra's algorithm can 
quickly calculate the Pareto set without the need to assign relative costs for 
the different modes. Then alternative cost functions can be evaluated on just 
the paths in the Pareto set or a post-optimal analysis conducted without ever 
having to rerun the algorithm on the network. 



4. Experimental Study 

In this section the weighted coloured-edge graph approach is applied to 
multimodal networks in two different scenarios. Firstly, the cardinality of 
J^uv is analyzed for complete graphs. Secondly, the approach is applied to 



10 



Table 1: Pareto set for network with 21 vertices and 51 edges. 
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Figure 3: Path 25 used for a post-optimal analysis 



a real multimodal network. This network corresponds to the transportation 
system of France and considers four transport choices. 

4.1. Number of processing paths and cardinality of J^uv 

The objective here is to identify general patterns for the number of pro- 
cessing paths and A4uv cardinality. In this test a weighted complete multi- 
graph is taken as input so that each analytical scenario is generated by fixing 
values for n = \V\ and k = \M\. Such a graph is characterized by hav- 
ing kn{n — 1) edges and the maximum number of possible paths |CP„t,| = 
(" ~ ^) /c-'"'"^j! for V u, which has factorial order O — 2)!). 

Specifically, the algorithm is run for complete multigraphs with k = 2, 3, 4, 5 
colours and n between 20 and 200 vertices. Random edge weights are gener- 
ated by means of a continuous uniform distribution of positive weights. 

Figure H] depicts the patterns followed by Aiuv cardinality. The figure 
uses a logarithmic scale for vertical as well as horizontal axes to demonstrate 
the average case polynomial behavior as n increases. Table [2] provides the 
numerical orders determined for different k values. These results demon- 
strate not only that the Pareto optimal set is calculated in polynomial time 
but that the resulting set requiring further analysis grows very slowly as a 
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Figure 4: Cardinality of A4uv for random weighted coloured-edge graphs with different 
number of colours k. 



fun ction of n. The results resemble ideas presented in iBentley et al.l (Il978h 
and lMiiller-Hannemann and Weihd ( 1200 ll ). suggesting the applicability of the 
model in real multimodal network scenarios, even when the networks are 
dense and without having to apply network reduction techniques or heuris- 
tics. 



4-2. Performance in a real multimodal network 

The approach is now tested on a large multimodal network. In this test, 
largeness is in the sense of number of vertices and edges. The selected net- 



Table 2: Order of processing paths and A4uv cardinality for several k values 



k 


Processing paths (p'su) 


Aiuv cardinality 
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0(^0.61) 
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work scenario corresponds to the multimodal transportation system of France 
being one of the largest networks in Europe. The multimodal network was 
obtained f r om v ector data information retrieved from a public GIS library, 



GeofabrikI fl2010h . 



The network dataset for each transport choice was firstly processed in 
ArcGIS to make it suitable for computation. ArcGIS is a Windows platform 
application for the analysis and processing of vector geographic information 
system data. This application has a network a nalysis e xtens ion that permits 
the identification of junctions and polylines (see iBurkd ( 120021 ) for a definition) 
in each transport system. In addition, ArcGIS also has a macro for the 
computation of the adjacency matrix for each system of junctions. 

Table [3] summarizes the number of junctions and edges for each transport 
mode as well as some statistics of the networks. All edge lengths are given 
in decimal geographic degrees. Four transport modes comprise the France 
transport system: road, rail, waterways and motorways. The road system 
mainly consists of primary roads. The rail system is comprised of common 
train lines disregarding subway and tram. Waterways are the channels and 
rivers used as transportation links. Finally, the motorway system of France 
includes toll roads and is considered a different mode of transport in its own 
right. As an illustration. Figure [5] depicts the France road system. 

The construction of the multimodal network requires assembling the data 
for the four network modes together. This task is accomplished by an ad- 
hoc algorithm coded in the Java language. The code basically takes two 
inputs. These are the adjacency matrix of each transport mode network 
and a list of minimum interjunction distances in each mode. The latter is 
built in ArcGIS by taking a mode junction dataset and applying the "join 
and relates" tool with respect to each other mode junction dataset. This 



Table 3: Characteristics of France Multimodal Network 



Modes Number of Number of Edge length 

Junctions polylines Maximum Mean Stnd. Dev. 

Roadways 53562 47660 0.868028 0.010674 0.032452 

Railways 18671 20083 1.280264 0.014966 0.046192 

Motorway 7720 7432 1.221951 0.033488 0.078485 

Waterways 17113 11635 3.238686 0.032070 0.095573 
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Figure 5: France roadway systeni. 



information facilitates the performance of a subsequent clustering procedure 
used inside of the ad-hoc algorithm. 

Two parameters need to be specified once the algorithm code is executed. 
A minimum clustering distance (this generates the network vertices) and a 
source vertex (a junction number). After entering this information, the Java 
code invokes the multimodal Dijkstra's algorithm, reporting at the end a list 
with the total number of optimal paths to each vertex together with two 
additional variables: the maximum number of paths found in a particular 
vertex and the average number of paths. 

This dataset was tested by assigning a cluster distances between 0.150 and 
0.115 decimal degrees (14 to 11 km). The resulting networks together with 
running times (minutes) and average number of optimal paths are shown in 
Table HJ The computations were performed on a standard desktop computer 
with dual core 2.93 GHz CPU and 8 GB of RAM that was set with the queue 
version of the multimodal Dijkstra algorithm. 

Although the networks shown in Table H] requiring longer runs of the 
multimodal Dijkstra's algorithm when the clustering distance is reduced, it 
cannot be disregarded that no constraints or reductions were required for 
obtaining the results in Table H] and the cardinality of the resulting Pareto 



15 



Table 4: Results of Prance Multimodal Network 



Cluster 


Number of 


Number of 


Running 


Avg. 


Max 




Vertices 


Edges 


Time 


Paths 


Paths 


0.150 


1501 


4216 


0.0182 


33.0930 


702 


0.140 


1869 


5218 


0.6953 


208.270 


3320 


0.130 


2343 


6280 


4.4010 


404.380 


14972 


0.120 


2948 


7696 


245.40 


2013.82 


37128 


0.115 


3108 


8018 


823.10 


3575.03 


46984 



sets are quite manageable for any further analysis. 
5. Conclusions 

In modelling multimodal networks, current approaches for determining 
shortest paths rely on applying application-specific constraints or heuristics 
to obtain tractable problems. This paper introduces a modelling approach 
that keeps the multimodal traits of a network by assigning discrete colour 
attributes to the edges, uses a partial order to obtain a Pareto set of paths of 
potential interest, and avoids the need for reduction techniques. Although a 
straightforward approach to modelling networks in which there are multiple 
transportation modes, it does appear to give a new perspective and truly 
general approach for multimodal networks. Another feature of this approach 
is that it results in a Pareto set, which can be further investigated without 
rerunning the algorithm. This opens the door to post-optimal analysis in 
MMN. 

The experimental study analyzing the Pareto set indicates that its cardi- 
nality is typically a low order polynomial for random uniformly distributed 
weighted coloured-edge graphs. Furthermore, the approach can deal with 
networks as large as the multimodal transportation system of Prance. 
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